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Abstract 

An elastic sheet lying on the surface of a liquid, if axially compressed, shows a transition 
from a smooth sinusoidal pattern to a well localized fold. This wrinkle-to-fold transition is a 
manifestation of a localized buckling. The symmetric and antisymmetric shapes of the fold 
have recently been described by Diamant & Witten (2011), who found two exact solutions 
of the nonlinear equilibrium equations. In this Note, we show that these solutions can be 
generalized to a continuous family of solutions, which yields non symmetric shapes of the fold. 
We prove that non symmetric solutions also describe the shape of a soft strip withdrawn from 
a liquid bath, a physical problem that allows to easily observe portions of non symmetric 
profiles. 

1 Introduction 

Localization of stress or strain is a phenomenon that appears in wide range of mechanical systems, 
like for example granular media 1 , fluid-saturated porous media^, twisted elastic rods^ or damage 
of materials'^. As extreme consequence, localization leads to the formation of singularities, like 
fractures^, break-up of liquid threads^ or d-cones in crumpled papers^, just to cite a few examples. 

In the field of elastic structures, the possible localization of the deformation of an axially 
compressed beam resting on an elastic foundation is a problem which has attracted a large interest. 
For this kind of systems, linear post-buckling regime is characterized by a sinusoidal shape with a 
well defined wavelength A, provided that the length of the beam is much larger that A™. However, 
it has be shown that the nonlinear post-buckling response exhibits a complex behavior, depending 
on the constitutive law of the foundation and thus on the type of bifurcation that occurs^,^. For 
instance, thin films on a soft elastic substrate display cascades of space-period doubling eventually 
leading to chaos^, as well as localized folding, delamination or creasing of the 

In the limit case of a thin sheet lying on a liquid surface, it has be shown that the sinusoidal wrin- 
kles appearing close to the buckling threshold (i.e for low confined configurations) spontaneously 
transform into a localized fold for highly confined configurations^. This phenomenon covers a 
wide range of length scales, spanning from lipid monolayers^ up to ice pack in the ocean^, and 
may lead to the delamination of the film from the liquid 18 . AudolyHSl has derived an amplitude 
equation for the envelope of the pattern, explaining the localization of the buckling and relating it 
to the one on a nonlinear elastic foundation. In a recent paper, Diamant & Witten^ have shown 
that the nonlinear equation governing the shape of the fold is integrable in the case of an infinite 
long sheet, and have proposed two exact solutions describing symmetric and antisymmetric shapes. 

In this Note, we show that symmetric and antisymmetric solutions are two particular cases 
included in a continuous family of solutions of the nonlinear equation. This family of solutions 
describes, in general, non symmetric shapes of the fold, all sharing the same total energy. While for 
the buckled floating sheet all the solutions are possible, we show that for a slight different system, 
consisting in the extraction of an elastic strip from a liquid bath, the exact shape of the strip is 
described by just one given solution among the family. 

2 A continuous family of solutions 

Consider the problem of a flexible and inextensible sheet of length 2L, width w and thickness e 
(with L ^> w 3> e), lying at the surface of a liquid (Fig. [TJ. The bending rigidity of the sheet is EI, 
E being Young modulus and I = we 3 /12 the second moment of the cross area of the sheet. The 
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Figure 1: 2D schematic view of a fold in a buckled floating beam. The angle 9 (between the tangent 
to the beam and the horizontal) as well as the horizontal x and vertical y positions, are expressed 
as a function of the arc- length s. The compressive dead load is noted P. 



density of the liquid is p, and g is the acceleration of gravity. We consider that the deformations of 
the sheet are confined in the (x, y) plane, which leads to study the equilibrium of a 2D elastic beam. 
The beam is parametrized by its arc-length s £ (-L, L). Horizontal x and vertical y positions of an 
infinitesimal element ds can be related to the angle 9 of the element with respect to the horizontal 
by: 



x'(s) — cos9(s) 
y'(s) = sin0(s) , 



(1) 
(2) 



where prime ' = d/ds denotes the derivative with respect to the arc-length. 

The linear response of this system to an axial compressive load P is a buckled sinusoidal 

/ a/4 

configuration, with a well defined wavelength A ~ ( ^ J , provided that L 3> A^. The linear 

response is valid only for a low confined beam (i.e. for small values of the displacement A at the 
ends of the beam). 

We introduce A as the reference length-scale, which is equivalent to take EI = pg = 1. We 
thus consider that all variables in the following are dimensionless. Here, as we are interested in the 
transition from a wrinkle state to a localized fold, we let L — ¥ oo and we make the assumption that 
all deformations vanish at infinity (y = 9 = 9' = at infinity). This leads to study the homoclinic 
orbit of the system, with the ends of the beam which represent the homoclinic connection in the 
phase space (0, 9')^. For these reasons, the fold appearing in the middle of the beam (which is 
flat otherwise) can be seen as a boundary layer solution. 

For an exhaustive derivation of the leading equations and a description of the nonlinear response, 
the reader can refer to the papers of Audoly^ or Diamant & Witten^. The total energy of the 
system, which is the sum of elastic, potential and compressive energy, writes: 



£(o,e',y,y') 



V( S ) 2 + ^(s) 2 costf( S ) 



P(l-cos0(s)) ds . 



(3) 



Functions 9(s) and y(s) are related by equation and the correct way to look for the equilib- 
rium equations is to find the stationary points of the lagrangian C = £ — J q(s)(y'(s) — sin#(s))ds, 
with q(s) a local lagrangian multiplier. This leads to the system of nonlinear equilibrium equations: 



9" / ( S ) + -9'(s) 3 + P9>( S )+y( S )=0 
y'(s) =sm9(s) 



(4) 



It is possible, by deriving the first equation in Q once, to write one differential equation for 
6(a): 



e""( s ) + l9'(s) 2 9"(s) + P9"{s) + sm9(s) = . 



(3) 
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Despite its complexity, Diamant & Witter^ were able to find two exact solutions to equation 
Q which satisfy the boundary conditions for 9 vanishing at infinity. These solutions are: 



s (s) — — 4arctan 
9 a (s) = — 4arctan 



csin(fc s) 
k cosh(c s) 

ccos(k s) 
k cosh(c s) 



(6) 
(7) 



with coefficients c and k linked to the load P by k = \J2 + P/2 and c = \/2 — P/2. Note that 
also — 9 S and —9 a are solutions. Here, 9 S is the solution which gives a symmetric shape and 9 a 
an antisymmetric shape of the fold. Indeed, 8 S is such that it vanishes at s — as well as all its 
even derivatives 6*", 6*"", ■■■ In contrast, 9 a is such that all the odd derivatives 6' a , 6'" , ... vanish at 
s = 0. 

However, these solutions are not unique, and it is possible to verify, by simple substitution into 
equation (|5|, that all the functions of the form: 



i(s) = — 4arctan 



csin(fc(s + <p)) 



k cosh(c s) 



(8) 



are equilibrium solutions of the problem (satisfying also the boundary conditions) . The parameter 
is a real number, and it appears only in the sine, so that it is not a mere change of variable 
for s. The fact that <f) can take any real values yields that equation ([8| is a continuous family of 
solutions. Hence, 9 S and 9 a are two particular solutions belonging to this family, recovered when 
= and 4> = n/2k, respectively. For <f> ^ mr/2k (n being an integer), the solution 0^ does not 
have any kind of symmetry, as we will show explicitly in the next paragraph. 

The shape (x^,^) of the fold can be obtained by integration of equations ([TJ-Q. Figure 
[2] shows, for two different loads Pi = 1.2 (left) and P2 = 0.5 (right), nine shapes belonging to 
the family of solutions. The corresponding values of (f> for each profile are given in the figure. By 
observing the shapes, one can easily remark that, beside the symmetric and antisymmetric profiles, 
there exists infinite non symmetric profiles. 

Diamant & Witten 20 already pointed out that symmetric and antisymmetric solutions have the 
same energy E . By substituting 9$ and into equation (13 L it results that £ — |V2 — P{2 — P) : 
the functional no more depends on <fi and all the profiles have the same energy. This means that 
all the solutions are neutrally stable with respect to each other. 



3 Selecting one non symmetric solution 

Among all the configurations parametrized by 4>, the question to know which ones are stable is 
still open. On the one hand, the images of Pocivavsek et alJ^ may suggest that some folds do not 
look completely symmetric, but the lack of experimental investigations on this topic does not allow 
for an exhaustive answer. On the other hand, it is possible to find a slightly different physical 
problem, in which the shape of the beam is exactly described by one given function 6^, the right 
value of cj) depending on a control parameter. 

In Rivetti & Antkowiak 2 ^, the authors study the equilibrium configuration of a 2D inextensible 
strip extracted from a liquid bath (Fig. |3j. The control parameter in the problem is the height of 
extraction, h, defined as the distance between the clamp and the level of the liquid. Parameter h 
can be tuned starting from zero (a completely flat strip) up to a maximal value at the moment of 
extraction. In this system, elasticity and hydrostatic pressure interact as for the buckled floating 
beam. However, in this problem there is no axial compressive load, but a pulling force at the end 
of the strip because of surface tension 7. In the case in which the length of the strip is L 3> A, the 
equilibrium configuration is given by the system of differential equations: 

9>"( S )+ 1 2 9> 3 -C 2 9'( S ) + y( S ) = 
y' \s) — sin 9(s) 
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Figure 2: Nine examples of solutions (x^, y,p) included into the family of solutions given by equation 
The solutions correspond to P = 1.2 (left) and P — 0.5 (right). The corresponding value of <j) 
is given for each profile. 




Figure 3: Schematic view of the extraction of a 2D elastic strip from a liquid bath. The strip is 
clamped at one end. The distance between the highest point of the strip and level of the liquid is 
noted h, and it is the control parameter of the extraction. 



with ( 2 = L 2 /X 2 representing the contribution of surface tension (L gc = Wj/pg is the capillary 
length)^. The minus sign for £ 2 , with respect to the compressive load P in equation is justified 
by the fact that surface tension pulls the strip. 

System of equations (j9| can be verified by all the solutions Q$ given by equation (JsJ) , if the 
substitution P = — C, 2 is done. Boundary conditions at s — > oo are identical to those of system 
Q and thus satisfied by 6$. However, the presence of a clamp requires to pay attention on the 
boundary conditions at s = 0. The condition of zero angle at the clamp, 9(0) = 0, can be settled 
by a simple change of variable s — > s — in equation ^ , which transposes 4> into the hyperbolic 
cosine: 

csin(fcs) 



-4 arctan 



k cosh(c (s — 4>)) 

Nevertheless, the control parameter explicitly appears in the last condition y(0) 



(10) 
h. Satisfying 
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Figure 4: Comparison between the experimental and the theoretical shape of an elastic strip 
withdrawn from a liquid bath. Here, the (x,y) theoretical profile (red curve) is obtained for 
h = 1.91A and ( = 0.24. leading to P — —0.06 and <fi = 1.33. The physical strip corresponds to 
a portion (the red dashed curve) of a non symmetric fold (the other portion of the fold is shown 
with a red continuos line). 

this condition leads to a selection of one specific value of <fi: 

y(0) = -j sin 00 (s) ds = 2^2 - P sech Q<^V2 - P^j = h (11) 

<j) = ^J=sech~ 1 [ — h | (12) 
y/2=P \2^/2^P) 

Hence, for any value of h, the shape of the elastic strip is given by one given solution 9$, which 
is a portion of a (generally) non symmetric fold. The Fig. [4] shows good agreement between the 
experimental profile and the corresponding solution obtained for P = — £ 2 = —0.06 and (f> = 1.33. 

From the analytical solution, it is therefore possible to obtain the value of the physical pa- 
rameters characterizing the extraction, as for instance the value of the curvature at the clamp, 
9'(0) — —h, or the vertical force needed to withdrawn the strip, F y = J °° y(s) cos 9(s) ds = 



4 Conclusion 

In this Note we considered the problem of a floating beam buckling under the action of a compres- 
sive load. In order to investigate the shape of the fold that spontaneously appears for high confined 
configurations, we presented a family of analytical solutions parametrized by a real parameter <j>. 
These new solutions predict non symmetric profiles of the fold, and they are a generalization of 
the symmetric and antisymmetric solutions already introduced^. We showed that all the profiles 
given by the solutions are neutrally stable with respect to each other, because they share the same 
total energy. However, the lack of an exhaustive study of the stability for all the solutions with 
respect to any perturbation does not allow to answer the question of which shape is selected when 
the fold forms. 

We proved that portions of non symmetric configurations can nevertheless be observed in a 
slightly different problem. Indeed, we related the problem of the buckling on a liquid surface to 
the problem of the extraction of a floating beam^. Although the equilibrium equations for both 
systems are very similar, boundary conditions in the latter are responsible for the selection of only 
one given solution. All the non symmetric profiles can therefore be observed by tuning the value 
of a control parameter appearing in the boundary conditions. 

One may ask whether the generalization introduced in this paper for the floating beam could 
apply to other soliton solutions, as in Sine-Gordon or Korteweg-de Vries equation. It is straight- 
forward to prove for both these cases that a solution keeps its validity if a ^shifting is applied like 
in this Note; however, as Sine-Gordon and Korteweg-de Vries equations are PDE, the (^-shifting 
merely leads to a change of variable for space or time, without perturbing the symmetry of the 
solution. Indeed, the intrinsic interest of the 0-shifting in the floating beam is that the equilibrium 
equation is an ODE, and 0-shifting effectively affects the symmetry of the solution. 
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This Note has proved that symmetric and antisymmetric solutions were not general enough. 
However, more general solutions than those introduced here could exist, because a direct resolution 
of the equilibrium equation has not be done yet. Moreover, a step forward in the comprehension 
of the problem may be the extension of the homoclinic solution to heteroclinic orbits in the phase 
diagram, i.e determine the exact shape of a buckled floating beam with a finite length. 
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